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Abstract 

Highlights of the results from ultrarelativistic heavy ion collisions at CERN-SPS 
are reviewed. In particular, I discuss how the experimental results indicate 
that a collective strongly interacting system has been produced, and what are 
the implications towards the Quark Gluon Plasma. The physical ideas behind 
measuring certain observables are introduced. The future program of high energy 
nuclear collisions at BNL-RHIC and CERN-LHC/ ALICE is also briefly discussed. 
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1. Introduction 

This talk is organized as follows: In the first, 
introductory section, I will briefly recapitulate the 
physics motivations for colliding heavy ions at 
ultrarelativistic energies. For the second section, I 
have selected some of the experimental highlights 
from ultrarelativistic heavy ion collisions (let me 
use an acronym URHIC) at CERN-SPS. For each 
topic, the basic idea will be explained first, followed 
then by a discussion of the interesting experimental 
observations. As URHIC at the SPS represent the 
highest energy nuclear collisions so far, they and 
their implications will have the main emphasis in 
this review. In the third section, I will briefly 
discuss the future program for heavy ions, BNL- 
RHIC and the ALICE experiment at CERN-LHC. 
In the fourth section, I will present my conclusions. 
For recent reviews on the same subject, see e.g. 



1.1. General facts 

The main goal of ultrarelativistic heavy ion 
collisions is to study the thermodynamics of 
strongly interacting matter. In colliding two heavy 
nuclei (A ~ 200) together at very high cms energies 
(a/s ~ 20...200...5500AGeV), the primary aim is 
to produce experimentally a new phase of matter, 
the Quark-Gluon Plasma (QGP), and - even more 
importantly - to observe the QCD phase transition 
from the QGP phase to the Hadron Gas (HG) 
phase. Cosmological motivation for URHIC is 
the fact that our Universe has undergone such a 
transition within its first microseconds. 



In general, although the basic motivation for 
such studies comes from the theory, it is fair to 
state that the field of URHIC is an experimentally 
driven one. In this field, it is very difficult to 
perform calculations from first principles, so - in 
addition to precision measurements - good, QCD- 
bascd phenomenology is needed. 

In relation to conventional high energy particle 
physics (HEPP), URHIC naturally share the goal 
of making high precision measurements. However, 
the two fields are orthogonal with respect to their 
physics goals: where HEPP aims for detection of 
new particles and new symmetries (Higgs, SuSy), 
URHIC reach for detection of a new phase of 
matter, made of the known QCD-quanta. In HEPP, 
simplicity is obtained through producing as few 
particles in the final state as possible, whereas 
in URHIC simplicity is obtained by producing as 
many particles in the final state as possible: only 
if the strongly interacting system becomes dense 
and large enough it can become collective and laws 
of thermodynamics become applicable. In other 
words, for URHIC the goal is to do elementary 
particle condensed matter physics. In this, and 
throughout my talk as well, the keywords are 
lifetime and volume of the produced, extended, 
system. 

1.2. QCD phase diagrams 

The theoretical foundation for searching the new 
QCD-phase of matter is given by the first principles 
calculations of lattice-QCD (for a review, see ||). 
From there we know that the QGP phase with 
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partonic degrees of freedom undoubtedly exists. An 
example of these calculations for the Equation of 
State (EoS, pressure as a function of energy density) 
of the QCD-matter is given in Fig. 0. The abrupt 
change in the energy density from the confined 
phase at low T to the deconfined phase at high T 
is a consequence of a first order phase transition in 
the pure gauge SU(3). With dynamical fermions 
but still with zero baryochemical potential /is = 
(which corresponds to zero net baryon number), 
the order of the phase transition has been seen 
to change with the number of quark flavours and 
with masses of the three lightest quarks. For the 
moment, inclusion of dynamical fermions can only 
be done in certain approximations and keeping 
fiB = 0. Calculations with /i^ / are not yet 
available, due to the complex valued actions. 





1 1 1 






/ / 




// / 


£ /T 4 


1: / 


3/4 s/J 3 




3p/T 4 






iJ 





Figure 1. Scaled energy density, entropy density and 
pressure vs. scaled temperature at the continuum limit 
m pure gauge SU(3) theory M. The horizontal 
dashed line is the ideal gas limit andthe vertical hatched 
band illustrates the latent heat. The figure is from 0. 

The critical deconfinement temperature T c at 
/i B = is about 265 MeV for the pure ST/ (3) 
glue, with quarks the estimates lie in the region 
T c = 140... 200 MeV §. While the deconfinement 
phase transition is studied through Polyakov loops, 
the phase transition related to chiral symmetry can 
be studied trough chiral condensates ($^ r ). In the 
lattice QCD calculations, chiral symmetry has been 
observed to be restored (i.e. (#\I>) — > 0) always at 
the same T c as the deconfinement transition occurs 

i- 

With URH1C, one explores the QCD phase 
diagram in the /is, T-plane as illustrated in Fig. |^. 
With increasing nuclear mass number A and with 
increasing cms-energy, more partons are liberated 
in the collision and the energy density of the initial 
QGP-system (immediately after the AA-collision) 
increases. With increasing initial energy densities 



and increasing initial temperatures, sketched by the 
arrows in Fig. 0, the lifetime and the volume of 
the plasma evidently increase as well, improving 
the chances for observing signals directly from the 
QGP. Baryon stopping has already been observed to 
decrease with the cms-energies J7|, and this trend 
is expected to continue at larger y/s. With less 
net baryon number stopped in the central rapidity 
region, which is the main region of interest for the 
QGP studies, also fiB will decrease, as illustrated in 
the figure. Ideally, it would of course be beneficial 
to get as close as possible to the theoretically best 
understood limit hb/T — > 0. Cosmologically, the 
interest lies also in this region: in the early Universe, 
the inverse of the specific entropy is tiny, ~ 10~ 9 . 

At very high cms-energies, the formed QGP 
can be said to be produced by "heating" : the 
QGP consists of gluons and quark-antiquark pairs, 
and the energy for creating them is provided by 
the high cms-energy. The QGP can exist also at 
T = at large values of \j,b- This would correspond 
to preparing the QGP by compressing the normal 
nuclear matter beyond the critical density, which is 
of the order of ~ 10 times the normal nuclear matter 
density 0.17 fm~ 3 . This may happen in the cores of 
neutron stars but as this region is not relevant for 
URHIC, I will not discuss it here. 
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Figure 2. An illustration of the conventional QCD 
phase diagram in the plane of chemical potential \lb 
and temperature T. The regions of the phase diagram 
the highest energy heavy ion collisions are probing arc 
sketched by the arrows. 

The current and future experimental program of 
highest energy nuclear collisions is also illustrated 
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in Fig. |2| In addition to the experiments I 
will discuss below, also an extensive intermediate 
and lower energy program for heavy ions (BNL- 
AGS, GSI-SIS, etc.) exists but in this talk I will 
concentrate only on the highest energy collisions. 
In the fixed target experiments at the Super 
Proton Synchrotron (SPS) of CERN, the highest 
beam energies have been Eb = 200v4GeV for 
32 S, and currently E b = 158AGeV for the Pb- 
beam in 208 Pb+ 208 Pb collisions. The corresponding 
nucleon-nucleon cms-energies are y/s — 20 GeV 
and 17.3 GeV. A lower energy run at the SPS 
with Eb = 40ylGeV will take place this year. 
In the near future, starting in December 1999, 
in the Relativistic Heavy Ion Collider (RHIC) at 
the Brookhaven National Laboratory (BNL), the 
collisions of highest energy and heaviest nuclei will 
be 197 Au+ 197 Au at Vs = 200,4GeV. In year 2005, 
in A Large Ion Collider Experiment (ALICE) at the 
Large Hadron Collider (LHC) at CERN, one will 
collide 208 Pb nuclei on 208 Pb at y/s = 5500 AGeV. 

As it is expected that baryon stopping nearly 
vanishes at the LHC energies, the experimental 
URHIC program will probe an extensive range of 
the conventional QCD phase diagram in the )Ib,T- 
plane. The critical question is whether the QCD 
phase transition can indeed be observed. Regarding 
the QCD phase diagram itself, there have been 
very interesting theoretical developments recently: 
at high densities, quarks may form Cooper pairs 
and a new, color superconducting, phase may exist 
[to). Most probably this happens at larger 
values oi hb than will be reachable experimentally 
by URHIC. On the other hand, the ^s,T-phase 
diagram may have a critical or tricritical point 
somewhere along the phase transition line in Fig. [|, 
depending on the order of the phase transition 
in full QCD - which depends on the masses of 
dynamical quarks and the effects of finite [Ib- By 
varying the cms energy it may be possible to trace 
this point down. For a recent review of these 
exciting theoretical developments, see |l(J. 

1.3. Space-time evolution: lifetime and volume 

The basic difference between URHIC and collisions 
of point-like (e + +e~) or small composite particles 
(P+P> P+p) is that an ideal system produced in 
URHIC is an extended and collective one with 
a large volume and a long lifetime. Once the 
system is produced, its space-time evolution cannot 
be controlled. In fact, the only experimentally 
controllable initial parameters are the mass 
numbers of the colliding nuclei and the collision 
energy. In addition, a handle on the impact 
parameter in each collision can be obtained by 



forming event classes of different multiplicities and 
transverse energies with a correlation to the energies 
observed in the zero-degree calorimeter. With these 
few controllable initial parameters, information of 
the whole space-time evolution of the system must 
be extracted from the various observables measured 
in the final state. 

The different stages in the space-time evolution 
of a strongly interacting system produced in an 
URHIC can be pictured as in Fig. pi At sufficiently 
high energies, the colliding nuclei are Lorentz- 
contracted thin disks surrounded by virtual clouds 
of partons. The original impact of the nuclear 
disks takes place almost instantaneously (within 
a transit time r tr = 2.R4./7) in a region around 
z,t ~ 0. Primary production of gluons, quarks 
and antiquarks occurs in the central rapidity region 
during some typical formation time To of the order 
of a fraction of fm/c, and r tr < tq. In a dense 
enough system the mean free paths of the QCD- 
quanta are much less than any typical homogeneity 
size V in the system: A^ fp <C V. A dense 
enough system reaches quickly a kinetic equilibrium 
with locally thermal quark and gluon momentum 
distributions. Chemical equilibration is preferable 
but not necessary to define a QGP, i.e. the density 
of quarks and antiquarks relative to that of gluons 
does not need to coincide with that in an ideal 
gas. The system will further evolve collectively as 
a QGP if it is interacting strongly enough and if 
the scattering times of quarks and gluons are much 
less than the inverse of the expansion rate (Hubble 
constant -1 ) of the system. For a recent discussion, 
see e.g. |J. 

In the most ideal case - in the limit of a 
very dense system - the evolution of the system 
after its formation can be described in terms 
of hydrodynamics. If the system were exactly 
longitudinally boost invariant, physics would be 
invariant along hyperbolas of constant longitudinal 
proper time r = \Jt 2 — z 2 [ jO] . In reality, however, 
already due to its edges the system is never perfectly 
boost invariant, which is illustrated in Fig. || by 
the curves bending backwards in time near the 
light cones. Notice also that depending on baryon 
stopping the temperatures and chemical potentials 
in the central rapidity region can be quite different 
from the ones at large rapidities. 

Analogously to the early Universe in the 
Big Bang cosmology, the produced QGP system 
expands and cools. Eventually the system ends 
up in a Hadron Gas phase. If the QCD-phasc 
transition is of first order, a stage of a mixed phase 
with both QGP and HG coexisting at a (local) 
critical temperature and chemical potential will 
appear in between. In the end of the HG phase, 
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Figure 3. Space-time evolution of a strongly 

interacting system produced in a very high energy A + A 
collision, projected in the plane of the longitudinal 
coordinate z and time t. Different stages of evolution 
and different types of observables are demonstrated. 



the average mean free paths of hadrons exceed the 
homogeneity size of the system, A^ fp > V, and 
the Hadron Gas decouples, i.e. freezes out. The 
hydrodynamical, locally thermal, picture ceases to 
apply at this point. In terms of hydrodynamics, the 
final state hadrons are emitted from the freeze-out 
surface. Notice, however, that the decoupling does 
not not take place instantaneously: as illustrated in 
Fig. ||, due to Lorentz-dilation effects, the regions 
near light cones freeze out typically later than the 
central regions. Furthermore, even within a certain 
spatial region decoupling is never an instantaneous 
process (as drawn in the figure for simplicity) but 
a dynamical one, determined by the expansion rate 
and the hadronic rescattering cross sections. 

The probes of strongly interacting matter can 
be divided into hadronic and electromagnetic ones. 
Typically, as illustrated in Fig. || the hadronic 
probes reflect the conditions at the freeze-out 
hypersurface, i.e. at the very last stage of the 
collective evolution. The electromagnetic probes, 
thermal lepton pairs and photons |l^, [H| I4| , in 
turn are emitted throughout the whole space-time 
history of the system; they decouple immediately 
due to their long mean free paths. The thermal 
electromagnetic probes have to compete with the 
primary production processes, e.g. Drell-Yan 
process for dileptons and prompt photon production 
for photons, which have the same origin in hard 



scatterings as in hadron-hadron collisions. Also 
the hadrons emitted from the freeze-out surface 
may decay electromagnetically and contribute 
strongly to the background for the electromagnetic 
probes. Because of the rapidly falling exponential 
momentum distributions in the thermal system, the 
momentum or mass scale of an electromagnetic 
probe is typically close to the local temperature 
at the time of emission. Therefore, the harder 
the probes are, the hotter and earlier system they 
typically probe. The probes coming directly from 
hard scatterings in the primary nuclear collision are 
the earliest ones and usually referred to as "hard 



probes" |15 . Some hard probes, like jet production 
(also in connection with direct photons) are strongly 
interacting but - thanks to the large momentum or 
mass scale Q ^> T involved - they are not absorbed 
in the QGP or HG. 
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Figure 4. Snapshots of an expanding system in 
longitudinal-transverse coordinate z, a;-plane, and in 
transverse coordinate x,y-plane. The black arrows 
indicate the longitudinal and transverse flow. The rest 
frame radius of a nucleus A is denoted by Ra- 

At very high energies, in a first approximation 
the produced system expands only longitudinally 
|p"lf . If the colliding nuclei were transversally 
infinite disks, this would be the case because 
no pressure gradients would be generated in the 
transverse direction. In reality, however, the 
colliding nuclei have rapidly falling nuclear density 
distributions ("edges"), and pressure gradients 
in the transverse direction are generated in the 
production and thermalization stage. Then, 
if the lifetime of the collective system is long 
enough, a strong transverse expansion develops |f6| . 
Consequently, at freeze-out the transverse size of 
the system can become clearly larger than the 
original nuclear radius Ra ~ 7 fm for A ~ 200. This 
is illustrated in Figs. ^. The quantitative details 
of evolution of the transverse flow depend on the 
actual densities and pressure gradients generated in 
the very beginning of the QGP phase, and on the 
actual EoS in the locally thermal system. The z, x- 
projection in Fig. [| demonstrates the longitudinal 
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and the transverse flow. Notice here that at 
fixed time t the different spatial regions in z 
at different stages of the evolution, as illustrated by Q 
Fig. | 

2. Highlights from the SPS 

In this section, I have picked up some of the 
highlights of the experimental results from URHIC 
at the SPS. I will concentrate on the observations 
that make URHIC very different from hadron- 
hadron collisions. In particular, I will discuss the 
experimental evidence of collective behaviour of 
strongly interacting matter. I will briefly review 
the anomalies observed in URHIC at the SPS: the 
strangeness enhancement from WA97, the low-mass 
dielectron enhancement from CERES/NA45, and 
the J/^-suppression from NA50, and discuss the 
physics implications of these experimental facts. 



158 A GeV/c Pb + Pb 



2.1. Transverse flow 

As discussed in the previous section, if the lifetime 
of the produced collective system is long enough, 
a strong transverse flow is generated. At freeze- 
out, the system emits particles from its decoupling 
surface according to the appropriate thermal 
distributions at each local freeze-out temperature 
Tf and baryochemical potential These 
distributions, however, must be folded with the 
(local) flow velocity in each emitting cell. Since 
the flow is a collective phenomenon, particles of 
different masses obtain the same velocity from 
it, and consequently the heavier the particle is 
the more transverse momentum it gains from 
the transverse flow. The collective transverse 
flow vt obviously therefore makes the transverse 
momentum spectra of heavier final state hadrons 
flatter than those of the lighter hadrons. In the 
mr-distributions dN/dm^, ~ exp[— tut /T s i ope ], the 
inverse slope parameter T s iopo can be expressed as 



1 



' slope 



Jfo + -m(v T ) 



(1) 



at the nonrelativistic domain p? <C m and as a 
"blueshift" formula 



^slopc T[ 




(2) 



at the relativistic region pr ^§> m 

In Fig. P jHj, also shown by Helstrup 
(WA97) in this conference, the observed m^-spectra 
of hadrons from NA44 |(J, NA49 @,|l],||, WA97 
pi, and WA98 pi, in Pb+Pb collisions at E b = 
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Figure 5. Mass dependence of the inverse slopes 
Tsiope of the mr-distributions. The fit is of the form 
dN/dm? ~ expf— mWTLi r, P J as compiled in pq| . The 
data is from 0^o|, 21 p3|,[24| and the figure is from 
the SPS-summary [!£]. 



158/lGeV, have been fitted to extract the inverse 
slope parameters T s \ ope as a function of the particle 
rest mass. Provided that the particles freeze out at 
approximately the same temperature, r s i op o should 
clearly increase with mass. This indeed seems to 
be the case, as seen in the figure, with a notable 
exception of fl, perhaps also that of 5. It has 
been suggested that due to their low interaction 
cross sections, these particles freeze out earlier and 
thereby benefit less of the transverse flow [^5| . 

The average transverse momentum has also 
been observed to grow with the size, as shown in 
Fig. ^. Obviously, while the general features are 
consistent with a collective picture, some details 
i.e. behaviour of fl (and S?), are not yet fully 
understood, and more work remains to be done on 
the theory. 

2.2. Asymmetric flow 

A feature in URHIC that in hadron-hadron 
collisions cannot in practice be addresses at all, 
are the non-central collisions, i.e. collisions with 
non-zero impact parameter b. For details of 
the asymmetric flow effects, I refer to a review 
p7| , but to get some intuition, let us look at 
Fig. |^, which is a snapshot of an A+A collision 
projected in the transverse (azimuthal) plane at the 
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Figure 6. The dependence of the inverse slope 
parameter on the size of the system as function of mass 
as observed by the NA44 Collaboration. The figure is 
from H. 



time of maximum overlap of the colliding nuclei. 
The generated pressure gradients will obviously be 
azimuthally asymmetric when b ^ 0. Transverse 
flow is generated by the pressure gradients in the 
transverse plane, so it will evidently be azimuthally 
asymmetric: (v x ) ^ (v y ). This in turn causes the 
azimuthal angle distributions of final state hadrons 
to be asymmetric. 
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Figure 7. A non-central collision projected in 
the transverse-coordinate plane, illustrating the initial 
azimuthal asymmetry in the forming pressure gradients 
(the arrows). 

To see the azimuthal asymmetry in the 
measured transverse momentum distributions, a 
Fourier analysis is performed, and a fit of the form 



F{4>) ~ 1 + 2vi cos (j> + 2v 2 cos 2</> - 



(3) 



can be used in each bin of rapidity. In practice, it 
suffices to consider the first two coefficients V\ and 
i>2- Figs. H illustrate the two types of asymmetry: 
purely directed flow with v\ ^ and v-i = and 
purely elliptic flow with vi ^ and V\ = 0. The 
measured angular distributions are then a mixture 
of these two effects. 

With the high precision data coming from the 
SPS, it has indeed now become possible to observe 
the few percent asymmetries in the azimuthal angle 
distributions. NA49 collaboration has measured the 
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Figure 8. Purely directed (at left) and purely elliptic 
(at right) flow patterns as the azimuthal asymmetries 
of the transverse momentum distributions. 



asymmetries in pion and proton distributions as 
functions of rapidity [fM|29[]. As an example, the 
recent results from [[29fl for the Fourier coefficients 
vi and vi for the pion spectrum are shown in Fig. [)| 
For the results from NA52, see Q and from WA98, 
see il. 
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Figure 9. The observed azimuthal asymmetries in the 
transverse momentum distributions of pions in Pb+Pb 
collisions at Ei = 158 ylGeV as function of rapidity and 
with impact parameters 6 < 11 fm, as measured and 
analysed by NA49 §§]. 

The existence of elliptic flow at the central 
rapidity (where the directed flow is zero) is a 
signature of pressure, and thereby collectivity, in 
the system. Once this signature is found, it is useful 
to scan the effect as a function of impact parameter, 
which controls the energy density reached. An 
abrupt behaviour as a function of b would be a 
direct signature of a QCD phase transition ]3^| . In 
the (preliminary) analysis p9| of the NA49 data no 
unexpected behaviour has been observed so far. 

2.3. Particle ratios and thermal models 

Above we saw how the experimental results support 
the existence of flow and collectivity in the system. 
The flow effects are expected to change the 
momentum distributions of final state hadrons but 
the total (full coverage 47r) yields of each hadron 
species should not, however, be sensitive to the 
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changes in the differential distributions. In fact, 
it can be shown that if all the fluid cells freeze 
out at the same temperature Tf Q and at the 
same baryochemical potential fi ! g everywhere, the 
cancellation of flow effects in the particle ratios is 
complete. Under these conditions, the underlying 
dynamics of the system can be ignored, and purely 
thermal models can be used to predict the particle 
ratios. The results should then depend only on 
Tf and fig, and by fitting the measured (full 
coverage) particle ratios with the predictions of 
thermal models, Tf Q and fJg can be extracted. 

Quite obviously, however, the particle spectra 
cannot be measured with a 4-7r-coverage but some 
extrapolation is always needed to get the actual 
47r-particle ratios. Secondly, the system is not 
likely to freeze out with the same Tf Q and fig 
everywhere, (in the forward and backward region 
the net baryon density is more than in the middle) , 
so when discussing a thermal model fit to the 
data, one effectively considers a substitution of 
the local quantities Tf a (x) and fi f g{x) by global 
averages (Tf Q ) and (fig). The good news, however, 
is that the hydrodynamical simulations |33[ show 
that only O(10 - 15%) deviations from (2&), (fi% 
are expected at the SPS energies. This implies 
that thermal model approach is a reasonable 
first approximation for extracting the freeze-out 
parameters Tfo and //g. For more discussion and 
references, see e.g. J^M. 
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Figure 10. Particle ratios as obtained in a thermal 
model Q with T fo = 170 MeV and fi% = 270 MeV as 
a best nt to the measured particle ratios from NA44, 
NA49, WA97, WA98 and NA50. 

As an example of a recent thermal model 
analysis, the results for particle ratios from are 
shown in Fig. [l^. The analysis is straightforward: 
a grand-canonical ensemble is used to describe 
a thermal system of hadrons at a certain global 



temperature Tf , in a volume V and at (global) 
chemical potentials fig, fis, fti 3 - Conservation 
of net baryon number, strangeness and charge, 
correspondingly, fixes V, fis and fii 3 , leaving Tf Q 
and fig as the parameters which are fitted to 
an extensive collection of data. In addition, the 
hadronic interactions must be modelled into the 
partition function, and one should also allow for 
decays as well (see for details). The best fit to 



the data from NA49 



WA98 



WA97 




NA44 |20j, 
NA~50 |4§, gives 



Tfo = 170 MeV and f7° = 270 MeV 



The results shown in Fig. |10| indicate that the 
measured particle ratios can indeed be quite well 
described by a system in thermal and chemical 
equilibrium. Thus a picture of a collective system 
in thermal and chemical equilibrium is emerging. 
The remaining differences between the model and 
the data are within the systematic uncertainties 
of the model and between different data sets 
I should also mention that there has been a lot 
of discussion whether the total strangeness is in 
chemical eq uilibrium or not | fl3| . In light of the 
analysis 1 35 , additional fugacities are not necessary. 



2.4- Identical particle inter •ferometry 

The main goal of identical particle interferometry, 
analogous to Hanbury Brown-Twiss (HBT) interfer- 
ometry of stellar objects H], is to extract a space- 
time picture of the system at freeze-out. By using 
the measured correlation functions in the momen- 
tum spectra of two (or more) identical pions, it is 
possible to estimate the emission volume (homo- 
geneity volume), lifetime of the system and dura- 
tion of the emission process (i.e. freeze-out). Also 
information on collectivity, i.e. transverse flow, 
is obtained. For more detailed reviews of HBT- 
theory and measurements in URHIC, I refer to 
here I will only describe the basic 

idea. 

Let us imagine an emission of two pions with 
momenta p\ and P2 where the sources are a distance 
R apart. Relative momentum is denoted by Q, as 
in Fig.[ll], left panel, and the total momentum by 
K. 

The two-particle correlation function is defined 

as 



C 2 (pi,P2) 



dN/d 3 pid 3 p 2 
dN/d 3 Pl dN/d 3 p 2 



(4) 



and it can be directly measured. The one-particle 
distributions (which are measured) are expressed in 
terms of space-time integrals of source functions 
S(x,p): dN/d 3 p! = J d A xS(x,p). The source 
functions contain all the information of the system 
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Figure 11. Left: Emission of two pions from two 
sources separated by R (a coordinate 4- vector). Right: 
The width of the corresponding (measured) two-pion 
correlation function C2(Q) reflects R: Q c ~ 1/R. 



at freeze-out, including the transverse flow. The 
non-trivial task is to unfold this information by 
using the measured correlation functions. 

The two-boson correlation function can be 
written in terms of the source functions as iff!! 



C 2 (Q,K) = 1 



\jd i xe^ x S{x,K)\ 2 
fd*xS(x, Pl )fd i yS(y,p 2 )' 



(5) 



The quantitative details naturally depend some- 
what on the physical assumptions made for the 
functional form of the source function but the ba- 
sic idea remains the same: the source radius can 
be determined from the width Q c of the measured 
correlation function, Q c ~ 1/i?, as illustrated in 
Fig. [0. 

In particular, it should be noted that unlike 
emission from the stars, which are static objects, 
emission from URHIC is a dynamic process. 
Consequently in HBT for URHIC, the spatial and 
temporal dimensions are non-trivially mixed, and 
several different radius parameters have to be 
defined. Usually K is chosen to be in the x, z 
plane with the beam in z-direction. Then the radius 
parameters are Rl, Rout and R S ide in the z-, x- and 
y-dircction, correspondingly. 

As an example of an HBT-measurement, let 
me show the one by NA44 in Fig. || @. The 
remarkable physics point here is that a typical scale 
0(20 MeV) ~ O(10 fm) is clearly seen in the high 
precision data of two-particle correlation functions. 
For the most recent experimental HBT analyses, see 
the results from NA44 §|, NA49 IDJ> WA98 & 
and the talk by Ganz (NA49) in this conference. 

As reviewed in | jf7| , when space-time informa- 
tion from HBT and complementary information 
from the my-distributions are combined, the freeze- 
out temperature of the system produced in Pb+Pb 
collisions is found to lie in the range 80 MeV < 
Tf < 120 MeV, a strong transverse flow is seen 
to develop, 0.35 < vt < 0.5, together with a 
large transverse extension R{ ~ 2i?pb, and lifetime 
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Figure 12. An example of the measured identical 
two-particle correlation functions in pion momentum 
distributions, measured by NA44 Js^] in Pb+Pb 
collisions. 



~ lOfm/c. In comparison with the thermal model 
results in the previous section, a clear hierarchy is 
observed: thermal freeze-out observed in the HBT 
and in the transverse momentum distributions takes 
place later, i.e. at lower T, than chemical freeze-out 
which is in turn reflected by the hadron abundan- 
cies. This is consistent with the expectation that 
the chemical reaction rates are much smaller than 
the elastic ones. 

2.5. Event-by-Event analysis 

Large multiplicities in single central events allow 
the study of fluctuations of different observables 
on an event-by-event (EbyE) basis. Also this is a 
feature of URHIC genuinely different from hadron- 
hadron collisions. Naturally, the primary goal of 
the EbyE physics is to observe diffcrcncies between 
the events with and without QGP-formation, as the 
QCD phase transition may cause large fluctuations 
which could show up in measurable quantities like 
multiplicities, (pt), particle ratios, etc in single 
events. Especially, if the QCD phase transition is of 
second order, one might directly observe the critical 
fluctuations associated with the phase transition. 
Also, if some more exotic physics takes place, such 
as formation of large enough domains of disoriented 
chiral condensates (DCC) j54|, the fluctuations 
from one event with DCC to another without could 
be seen. 

Fig. |l3| shows the recent EbyE results from 
NA49 [pSf ! So far no evidence of non-trivial 
dynamical fluctuations, i.e. fluctuations beyond 
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Figure 13. Event-by-event fluctuations of single event 
(Pt) (left) and single event K/-7r-ratio (right), measured 
by NA49 in Pb+Pb collisions at the SPS. The dynamical 
event-by-event fluctuations are < 1% in (pr) and < 
4.9% in K/tt. The figure is from EE). 



finite statistics and variations of impact parameter, 
has been observed in the Pb+Pb collisions at 
SPS: the dynamical event-by-event fluctuations are 
less than 1% in single event (px) and less than 
4.9% in K/7r-ratio (at 90 % confidence level). It 
will indeed be very interesting to see what the 
corresponding results will be in URHIC at RHIC 
and LHC/ALICE. 

2.6. Strangeness enhancement 

The observables discussed in the previous subsec- 
tions support the interpretation of the the system 
produced in Pb+Pb collisions at the SPS as a col- 
lective strongly interacting system. Let me now 
move on to other experimental facts, namely the 
observed anomalies. First, let us have a look at 
the strangeness enhancements observed. For more 
detailed reviews and for more complete list of refer- 
ences to the measurements of strangeness produc- 
tion in URHIC, see e.g. ||||. 

In comparison to hadron-hadron and e + +e _ 
collisions, a global enhancement of strangeness 
production has been observed in URHIC. At the 
SPS, the first observation of this was made by 
NA35 @ in S+S collisions at E b = 200AGeV. 
The global enhancement can be quantified by 
defining a parameter which counts the strange 
quark-antiquark pairs produced relative to the uu 
aid produced B9[ 



A 



AA 



2(8 + s) 



(u + u) + (d + d) 



2\P p 



(6) 



This indicator is plotted in Fig. [IJ for S+S 
and S+Ag at E b = 200AGeV based on the 
data from NA35 fl60|, and for Pb+Pb collisions at 



Eb = 158ylGeV based on the measurements by 
NA49 |3(], ^3l| and compared with the corresponding 
compilation of data from e + +e _ , p+p and p+p 
collisions (for the refs, see |5!|). The observation 
is that \f A rj 2Af p,ee : the global strangeness 
enhancement is thus a unique feature of URHIC. 



,0.6 



0.5 



0.4 



0.3 



0.2 



0.1 



o pp collisions 

A pp collisions 

□ e + e~ collisions 

• A+B collisions 



LO 



I 10" 
V s (GeV) 



Figure 14. The global strangeness is observed to be 
enhanced by a factor ~ 2 in A+B (S+S, S+Ag and 
Pb+Pb) relative to e + +e~, p+p and p+p collisions. 
The figure is from ]Hs| ]. 

In addition to the global enhancement, 
also specific enhancements in the production of 
K, K, A, A, H, S, il, O and 4> at midrapidity in 
Pb+Pb collisions relative to that in p+Bc collisions 
have been reported by the WA97 |||,||], NA49 
||,|| and NA50 || collaborations. The multi- 
strange hadron yields, in particular, are strongly 
enhanced in central Pb+Pb collisions, as seen in 
Fig. |l5], where the yields per participants relative to 
pBe, measured by WA97, are plotted as a function 
of the average number of par ticip ants in p+Pb and 
Pb+Pb collisions at the SPS ]62f. For more details, 
see the talk by Helstrup (WA97) in this conference. 

Notice especially that the number of negative 
hadrons scales with the number of participants, and 
relative to the charged hadrons the enhancement 
in £1~ + il + yield is roughly a constant factor 10. 
With the 1998 run, the Pb+Pb event statistics 
was practically doubled, and WA97 will soon cover 
the region below 100 participants (see Helstrup 
and J62[). It will certainly be very exciting to 
see whether there is a threshold in the number 
of participants (centrality) after which the strong 
specific strangeness enhancement takes place. 



10 



+ 
a. 
o 



MO 



B 
a 
a. 



-r-i--|-} h" 

— 


1 lJ- 

£2 +£2 


! 

pBe pPb PbPb 


~tl : 

pBe pPb PbPb 



250MeV<M, 



; 700 MeV, usually referred 



10 



10 



101 



10 



10 



10 

<N na „> 



Figure 15. The specific strangeness enhancement as 
measured by WA97. The yields per participant per unit 
rapidity relative to p+Be as a function of the number 
of participants (centrality) (62) . 



What are the implications of these experimen- 
tally observed facts? If the produced dense system 
consists of uncorrelated quarks (and antiquarks), it 
is intuitively clear that the more a hadron con- 
tains strange quarks (or antiquarks) the more it 
benefits of the global strangeness enhancement in 
hadronization. This picture is quantified in statisti- 
cal hadronization (quark coalescence) models J66[ , 
which seem to work for A+A collisions but not for 
p+A collisions. On the other hand, hadronic ki- 
netic models do not reproduce the data on multi- 
strange hadrons due to high mass-thresholds and 
low chemical reaction cross sections. Also ap- 
proaches with string breaking followed by hadron 
rescattering models (see |52| for the actual compar- 
ison and for further references) cannot reproduce 
the centrality pattern observed. The implication of 
the observed strangeness enhancement therefore is 
that the quark degrees of freedom are indeed essen- 
tial, and that strangeness must have been produced 
very early on in the collision, in the pre-hadronic 
stage. 

From the thermal models (Fig. [To]), we 
remember that the particle ratios, including those 
with multi-strange hyperons are close to chemical 
equilibrium values. This, together with the 
conclusion of the importance of the quark degrees 
of freedom, points towards thermalized QGP. 

2.7. Low-mass e + e" enhancement 

Another very interesting anomaly observed 
in URHIC at the SPS is the excess of 
electron-positron pairs in the mass region 



to as the "low-mass" region. The excess in e + e 
production was first observed by NA45 /CERES 
in S+Au collisions at E b = 200AGeV @,|6|0, 
and a similar excess was observed in by 
the HELIOS-3 collaboration in S+W collisions at 
E b = 200AGeV Q, and by NA38 @ in S+U 
at E b — 200 AGeV. The low-mass excess has 
been confirmed by CERES for heavier systems in 
Pb+Au collisions at E b = 158 AGeV @j7§. 

The reference data is given by the e + e~ 
production in p+Be and p+Au collisions at E p = 
450 GeV: no enhancement is observed there, the 
data can be explained by including contributions 
from all the known hadron decays. The mass 
distributions of e + e~-pairs measured in p+Au 
by CERES Q and scaled by the measured 
charged particle multiplicity (within the CERES 
acceptance) are shown in Fig. n6| The cocktail of 
the several different sources of €^e~ -pairs is shown, 
and it is seen to reproduce the measured data within 
the estimated errors (the shaded band) quite well. 
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Figure 16. Inclusive mass distribution of e + e~-pairs in 
p+Au collisions at 450 GeV scaled with charged particle 
multiplicity |7^] . The curves are the known sources from 
hadron decays. 

However, when one extrapolates the cocktail- 
plot to S+Au and Pb+Au collisions, a clear 
enhancement is observed, as shown in Fig. 
for the mass distributions of e + e~-pairs scaled by 
the charged particle multiplicity. In the figure, data 
sets from '95 and '96 are compared with the sum 
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of the expected contributions from hadron decays. 
In the cocktail-plot the particle ratios are taken 
from a thermal model fitted to measured ratios in 
Pb+Pb (T fo = 175 MeV and p B = 270 MeV, figQ, 
and the ?y-and pr-distributions follow the measured 
systematics in Pb+Pb collisions. The enhancement 



in the mass-region 250 MeV < M e 



700 MeV is 



observed to be 2.6 ± 0.5(stat) ± 0.6(syst) for the 
'96 data set, and 3.9 ± 0.9(stat) ± 0.9(syst) for the 
'95 data set j73|]. The excess has been observed to 
concentrate atlow pair transverse momentum. 
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Figure 17. Low-mass enhancement of e + e 
Pb+Au collisions at E b = 158,4GeV [fFl. 
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Interestingly, the excess scales more strongly 
than linearly with charged particle multiplicity, 
suggesting that it is due to a medium-effect. The 
shape of the mass-spectrum is not consistent with 
7r + 7T~ annihilation in free space, either, which again 
supports a medium-effect as the origin of the low- 
mass enhancement. There has been considerable 
theoretical activity in trying to explain the effect. 
For recent reviews, see H, [75|[72|]. The suggestions 
vary from a collisional broadening of the p-meson 
(i.e. p has a shorter lifetime in a dense hadronic 
medium, which leads to a broader peak in the mass 
distribution) and in-mcdium modifications of the p 
spectral function |l4|, [75| to a possible change in the 
mass of p ]76| . 

Let me also show the dielectron results from |77[ 
based on a hydrodynamical simulation where the 
whole space-time evolution of the system is taken 



into account [[78). In this approach, the total yield 
of dielectrons (and dileptons in general) consists of 
the pairs from the decays of resonances originating 
from the freeze-out surface (see Fig. ||) , and of the 
thermal dielectrons emitted from the dynamically 
evolving fireball throughout its whole space-time 
evolution. The measured hadronic momentum 
spectra are first reproduced to constrain the hot 
initial state and the EoS simultaneously j78|. After 
this, the emission of thermal e + e _ -pairs can be 
predicted by using the thermal rates 79, 7^ at 
each local T and ps- The e + e~ contribution from 
the resonance decays alone is shown in Fig. |l8| by 
the upper panel. Adding the thermal pairs gives 
the lower panel. We observe that that thermal 
emission dominates the total yield of e + e~-pairs in 
the region of the experimentally observed excess. 
Let me also remind that the low-mass enhancement 
is not considered to be a direct signal of QGP: 
contribution of thermal emission from the QGP- 
phase is negligible in the region of the low-mass 
enhancement (see |77j). 

Although theoretically still not fully resolved, 
it is evident that with the low-mass enhancement 
of e + e~-pairs measured by CERES, we are having 
a first glimpse into the dynamics of mesons in 
baryon-rich environment exceeding a temperature 
T ~ 120 MeV. In the near future, CERES will 
collect more data, also with the lower energy beams 
(E b = 40AGeV) and include an additional TPC 
in the experimental setup ]7q| . With improved 
statistics, and a better mass-resolution and signal- 
to-background ratio, CERES will certainly shed 
more light to the details of low-mass dielectron 
excess as a signature of a hot hadronic gas. 

Regarding the observed dilepton mass spectrum 
in general, let me note that also the intermediate 
mass dimuons, 1.5 GeV < < 2.5 GeV, show an 
excess relative to the conventionally known sources, 
as observed by NA38 in S+U fof, and HELIOS- 
3 in S+W collisions at E b = 200AGcV M, and 
NA50 in Pb+Pb at E b = 158 AGeV IsTl . It will be 



interesting to see whether this excess could also be 
explained by production of thermal dileptons Q . 

2.8. J/ty suppression 

Of the anomalies observed in URHIC at the SPS, 
the one that has generated most excitement, is the 
anomalous suppression of J/ty in central Pb+Pb 
collisions at E b — 158 AGeV, observed by the 
NA50 collaboration g, Q 11 II ■ The idea of 
studying the suppression ofJ/4' production as a 
signal of the QGP originates from Satz and Matsui 
in 1984 pM]: formation of J/ty bound state should 
be very efficiently Debye screened in QGP. For 
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Figure 18. Mass spectrum of 
with charged particle multiplicity 
a hydrodynamic simulation |W|. 



e + e~-pairs scaled 
as predicted by 
Upper panel: 

contribution from the hadron decays. Lower panel: 
total yield, which includes the thermal emission. 



details of the NA50 measurements and analysis, 
see |87| and the talk by De Falco (NA50) in this 
conference. For a more detailed overview of the 
theory of J/\l/ suppression, see e.g. 90,[)l| and 
the talk by Nardi in this conference. 

The key question naturally is relative to what 
the J/4" suppression appears in the heaviest URHIC 
systems. As shown in Fig. |l^ H|, a clear and 
smooth suppression pattern is observed from p+p 
and p+d collisions to p+A and S+U collisions, as 
measured by the NA38 and NA51 collaborations 
|p2f : B flfl a( J/ 1 ^) is decreasing as a function of A*B. 
For the comparison, the measurements have been 
rescaled to the same \/s. A scaling B^^a(J ~ 
(AB) - 92 is found |84| . As seen in the figure, 
the J/ty in Pb+Pb collisions clearly deviate from 
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Figure 19. Suppression of J/\& production as a 
function of AB: normal nuclear suppression from p+p 
to p+A and to S+U, and anomalous suppression in 
Pb+Pb. The figure is from NA50, ref. Q. 

At the same time, the Drell-Yan (DY) dilepton 

§duction at is not observed to be suppressed 
| a(DY) cxp /(7(DY) th - AB. Obviously then 
the suppression of J/^f comes from the final state 
interactions, i.e. from absorption of the produced 
bound state within normal nuclear matter (the 
beam and target). The nuclear absorption, often 
called as "normal" suppression^ can be understood 
in Glauber-type models, see 



96 for details. For 



discussion of the absorption by hadronic comovers, 
see [P5[. The observed nuclear absorption pattern 
(Fig. |19|) can be reproduced by using absorption 
cross-sections of the J/M/ (or its pre- hadronic state) 
of the order of 6... 7 mb p4, p6| 

NA50 has also carefully studied the suppression 
as function of transverse energy Et- In each Pb+Pb 
collision, the produced transverse energy and the 
energy observed in the zero-degree calorimeter 
correspond to an impact parameter b. The 
produced cc-pair bound state which is to form a 
J/\E' first has to go through the rest of the projectile 
and the target, so for each b (Et) there is an average 

| The Drell-Yan computation here assumes that any nuclear 
effects in the parton distributions | p^ ] are negligible within 
the errorbars. 
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nuclear path length L which the pre-resonance has 
to survive. The measured Drell- Yan pairs serve as 
a non-suppressed background in each £V-bin, i.e. 
at each 6, and B w <r(J/$)/(r(i)y) can be plotted 
both as a function of Et and as a function of L, 
as shown in Figs. M, |l] and || f|. The DY 
cross section <7(£>Y)7against which _B MM (7(J r /\E') is 
compared, is computed at 2.9 GeV < M < 4.5 GeV 
by using the measured DY-pairs at M > 4.5 GeV 
to get the correct normalization (Note that equally 
well one could directly use the measured Drell- Yan 
pairs; this has also been done by NA50 but is not 
shown here). 
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Figure 20. The suppression of J/\l/ over Drell- Yan 
pairs in 2.9 GeV < M < 4.5 GeV as a function of the 
average nuclear path length L of the cc pre-resonance 
[§6|, for NA38, NA50 and NA51 data. The NA50 
Pb+Pb data (filled circles) is from the 1996 run. 



In Fig. 20 



the filled circles at largest values 
of L represent the most central Pb+Pb collisions, 
and the ones at smallest L correspond to impact 
parameters (b) — 10.8 fm |36|. The peripheral 
Pb+Pb are seen to coincide with lighter systems, as 
one would expect to. It is also interesting to divide 
out the normal nuclear suppression and replot the 
suppression, as done in |86| , from where I have 
borrowed also Fig. 21. Notice here the importance 



of the p+p and p+d points (not shown in this 
figure) in determining the amount of normal nuclear 
suppression of J /if? at large values of L. The filled 
circles show the results from the new minimum bias 



analysis of NA50 p6[ , as discussed by De Falco 
at this conference. The abrupt drop seen in the 
figure for the J /if? suppression in Pb+Pb is the 
best candidate so far for a threshold behaviour in 
URHIC at the SPS. Naturally, the hope is that this 
would be the first direct observation of the QGP 
formation in the central Pb+Pb collisions. 
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Figure 21. The same as in the previous figure but 
the with the normal nuclear suppression divided out 
H. The Pb+Pb data is from 1996 run, and also the 
new minimum bias analysis is shown (the filled squares). 
In comparison with the previous figure, notice here the 
different scale in L. 

Fig. from [g6| shows the suppression of 
J /if? over Drell- Yan directly as a function of the 
measured Et, as obtained in the two different 
analyses of the data from 1996 run. The 15 Et bins 
correspond to the 15 bins in L in the previous figures 
and the solid line represents the expectation for the 
normal nuclear suppression as shown in the previous 
figures. Here, one should again keep in mind the 
importance of the lighter systems in determining 
the solid line. The suppression pattern is, naturally, 
the same as in the previous figures. 

Interestingly, as recently reported by NA50 |3?J , 
there may actually be another drop in the Jpf? 
suppression pattern at the largest values of Et- 
The latest (preliminary) analysis of the 1998 data is 
shown in Fig. |2^ from ref. [^7| . As discussed in [g7| , 
the second drop was not observed in the data from 
1996 run, where a thicker target was used, because 
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Figure 22. The ratio vj/v/vdy at 2.9 GeV < M< 
4.5 GeV in Pb+Pb collisions as a function of Et 
The 1996 data is analysed in two different ways, see |8(J 
for details. 



of a bias from re-interaction events. This problem 
is removed in the 1998 set up. The possibility for 
a second drop is indeed very exciting, and it would 
support the prediction in [ p7| of melting of different 
cc bound states at different energy densities. 

For a future improvement, the possibility to 
add a new vertex detector into the NA50 setup is 
discussed. This suggestion, called NA6z, was also 
presented in this meeting by Shahoyan. The benefit 
of this would be that open charm contribution 
could be measured more accurately, and this would 
serve as a very useful background both for the J/^f 
suppression measurements and for the intermediate 
mass muon pairs. 

For a compact review of the theoretical 
interpretations of the observed J/^ suppression, 
let me suggest [B9), see also the talk by Nardi in 
this conference. As one is discussing formation and 
interactions of a QCD bound state, perturbative 
calculations are hard to do. Consequently, different 
models have been presented, see e.g. [pl| p5j, |9§|, 

sTlikc the 



100, 101, 901 M. Several questions, 



energy depend ence of t he dissociation cross section 
of the cc-pair [^0l|l0l , and applicability of QCD 
factorization in the bulk production of J/'f' in 
AA-collisions have been and are still discussed. 
Even if all the details are theoretically not yet 
fully known, there seems to be a general consensus 
on that while the normal nuclear suppression of 
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Figure 23. Preliminary analysis of the 1998 data for 
suppression of J/\P and comparison with the analyses 
of the 1996 data. The figure is from @. 



J/\P can be accounted for by a more conventional 
hadronic picture, to reproduce the anomalous 
suppression, especially with the observed abrupt 
behaviour - possibly even with two drops - an 
additional ingredient very efficient in destroying the 
cc bound state is needed. Putting this together 
with the implications of a collective system and of 
the importance of partonic degrees of freedom, as 
discussed in the previous sections, it is very likely 
that we are now finally looking at a deconfined 
phase of QCD-matter. 

3. The future: RHIC & LHC/ ALICE 

3.1. RHIC: Relativistic Heavy Ion Collider 

Until now the fixed target experiments at the 
SPS at V« ~ 20AGeV have been the URHIC of 
highest collision energy. The next step will be 
taken by the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory, where 197 Au- 
beams will be collided onto 197 Au-beams at y/s — 
200/lGeV, and luminosity L ~ 10 26 " 27 cm -5 s -1 . 
RHIC will also collide proton-beams with proton- 
beams at y/s = 500 GeV and L = 10 31 cm^s -1 , 
and make p+ A and A+B collisions with various A 
and B. Collisions with polarized protons will also 
be done. At the time of the this meeting, beam tests 
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were being done in the two rings, with the aim to g et 
the collider in operation in December 1 998(§| [ 104 1 . 

Without going to any details, let me here just 
list the heavy-ion experiments at RHIC, together 
with their main ph ysic s goals. For a more complete 
review, see e.g. [105 . The main goal of RHIC 



will be to study the strongly interacting matter 
produced at URHIC, and to observe the QGP 
through various independent signals. Thanks to the 
varying beams and cms-energies, RHIC will provide 
us with an extensive y/s- and AB-systematics, 
which will be important in comparing the events 
with and without QGP formation. 

There are four major heavy-ion experiments 
at RHIC, involving 0(900) scientists and O(80) 
institutions. There are two "large" ex perim ents, 
STAR (Solenoidal Tracker at RHIC) (l06) and 
PHENIX (Pioneerin g Hig h Energy Nuclear Inter- 
action Experiment) 107 |, and two "small" ones, 
BRAHMS (Broad Range Hadron Magnetic Spec- 
trometer) |0§ and PHOBOS |T09 |. 

The main elements of the STAR experiment 
are silicon vertex tracker system (SVT), time 
projection chamber (TPC), forward radial-drift 
TPC (FTPC), time of flight measurement system 
(TOF) and electromagnetic calorimeter (EMC). 
The primary physics goal of STAR is to study 
strongly interacting matter at high densities in 
central A+B collisions. This will be achieved by 
measuring hadrons over a large solid angle, and 
measuring the global observables on an event-by- 
event basis. STAR will also study peripheral A+B 
collisions and polarized p+p collisions. 

The main building blocks of the PHENIX 
experiment are multiplicity vertex detector (MVD), 
ring-image Cherenkov counter (RICH), time- 
expansion chamber (TEC), TOF and EMC. 
PHENIX will mainly concentrate on measuring 
leptons, photons and hadrons to search for the QGP 
signatures. 

BRAHMS is a forward and mid-rapidity hadron 
spectrometer. The physics goal is to make a 
systematic study of particle production in AA 
collisions from the peripheral to the most central 
impact parameters. BRAHMS should basically 
answer the question of where the baryons go. 
PHOBOS is a compact multiparticle spectrometer 
which will measure single particle spectra and 
correlations between particles with low transverse 
momenta and to characterize events using a 
multiplicity detector. 



Some delays can now be foreseen. 



3.2. ALICE: A Large Ion Collider Experiment 

The ultimate step in the cms-energy of URHIC 
will be t aken by A Large Ion Collider Experiment, 
ALICE JIhJ at the Large Hadron Collider (LHC) 
at CERN in 2005. ALICE will be the only detector 
at the LHC dedicated for URHIC. 

The heaviest beams colliding at ALICE will be 



pb+ 208p b at ^ _ 5 5 ^TeV. The luminosity ii 



208j 

these collisions, L — 2 • 10 27 cm _2 s _1 , corresponds 
to some 10 4 interactions/s. ALICE will also 
perform p+p and p+A collisions to collect reference 
data for Pb+Pb, and it will also make A+B 
and A+A collisions with lighter nuclei to vary 
the produced densities of the QGP. According to 
the present running scenario, the running time for 
heavy ions will be about 10% of the LHC year. 

In comparison with SPS and RHIC, there are 
definite physics advantages in increasing the cms- 
energy as much as possible: 

• net baryon density near y ~ is smaller; 
[Lb -C T, and the EoS approaches that of the 
theoretically better understood fis = case 

• initial QGP energy density is higher 

• lifetime of the QGP phase and of the whole 
system is longer 

• freeze-out volume is larger (in a comoving 
sense, the system never decouples simultane- 
ously everywhere) 

• one moves into the applicability region of 
perturbative QCD (pQCD): First, there will 
in general be higher rates for harder QCD- 
probes of the QGP, such as high-pT jets, direct 
photons at high px, large mass dileptons, 
and upsilon production, and, second, the 
initial energy densities can be more reliably 
estimated based on perturbativ ely calcu lable 

see 



quark and gluon i/p-production [111, 112 



also the talk by Hammon at this conference, 
and, the approach of refs. [ 1 14 1 . Let me 
point out that the pQCD calculation has now 
been do ne c onsistently in next-to-leading-order 

pQCD Jn|. 

In this conference, the detector overview and 
physics capabilities of ALICE were reviewed by 
Morsch (ALICE). The main goal of ALICE is 
to study the behaviour of strongly interacting 
matter through observing various independent 
QGP-signatures simultaneously. A recent revi ew o f 
the physics at ALICE can be found in e.g. |115|. 
Let me, however, list some examples below: 

• The global observables, such as charged parti- 
cle multiplicity at central rapidities, dN c h/dy, 
and global transverse energy distributions 



dEr/dy (da/dEx) will be affected by the ini- 
tially produced densities (initial transverse en- 
ergy is not an observable unlike the final Et) 
and subsequent pdV work done by the collec- 
tive QGP-HG system against expansion [112]. 
Event-by-event fluctuations will be measured 
in many observables, such as forward-energy, 
charge particle multiplicity, HBT-correlations, 
single event (pr), particle ratios, N c h/N 1 . In 
this way ALICE is expected to find the possible 
exotic events. 

The transverse momentum spectra 
and particle abundancies will be mea- 
sured for all possible particle species 
(ir, 77, u, (j),p, K, A, 3, fl, D, d, t, a, ...). These 
will reflect the flow effects and the freeze-out 
as discussed previously. 

Identical particle HBT-correlations will be 
measured to give information of the dynamical 
evolution, size and lifetime and phase-time 
structure of the freeze-out process. 
Open charm and open beauty measurements 
will provide a normalization for J/ty- and T- 
suppression. As hard probes, cc and bb are 
produced very early on in the AA collision. 
Consequently, they will be very sensitive to the 
early and hot QCD phase of matter 
Energy-loss of high-py partons in QGP 
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can be studied via jet (and leading particle) 
measurements at ALICE. In this conference, 
the recent results on finite systems were 



discussed by Levai [117]. 
Muons at 2 < rj^ < 4 will be measured and 
will be studied. Especially interesting 
will be the dimuons coming from the decays 
of J/ty and T (if there are any J/W,T to 
decay!), as these will be directly probing 
the deconfinement. Also intermediate mass 
muons will be very interesting, due to the 
contributions from several sources, such as 
the Drell-Yan process, decays of D - and B- 
mesons and thermal emission [119]. If e + e~ 
pairs will also be identified, as lately suggested, 
modifications of the p, u>, 0-meson properties 
related to the low-mass behaviour of e + e~ can 
be studied in detail. 

Hard photons will also be an important tool 
to study the QGP. Especially, in order to see 
any window for thermally radiated photons, 
the background of direct photons and decay 
photons needs to b e ful ly understood. As 
discussed for RHIC ]118| , the direct photon- 
jet events can possibly also be used to study 
energy loss of partons in QGP. In these 
calculations, and in calculations of any hard 
probes of QGP, information of nuclear parton 



distributions and their scale evolution |93| arc 
needed. The data from p+^4 collisions will 
also be used to constrain the nuclear parton 
distributions. 

To compare some global characteristics of 
URHIC at SPS, RHIC and LHC/ ALICE, I have 
prepared the table 1. I would especially like to 
draw your attention to the expected long lifetime 
of the QGP-phase at the LHC which is of the same 
order of magnitude as that of the whole system 
at the SPS. Also, even though the system will 
not freeze out simultaneously everywhere in space 
(roughly simultaneous decoupling may happen in 
the transverse direction, except for the very edges, 
but longitudinally there will be strong time dilation 
effects at the LHC), I have defined a comoving 
volume by counting roughly how many units of 
space-time rapidity are available along a proper 
time curve. Notice how dramatic the increase in 
the "size" is when going from the SPS to the 
LHC: the rest frame volume of a lead nucleus is 
Vp b = 1150 fm 3 , so 400Vp b ~ 5 • 10 5 fm 3 , and 
compared with the nuclear size parameters, one is 
talking about almost "macroscopic" quantities of 
strongly interacting matter at the LHC! 

To conclude the LHC-section for URHIC, let me 
also mention that the possibility of measuring J/^, 
T and Zq production in Pb+Pb collisions in the 
CMS det ector at the LHC is being discussed, see 
e.g. p4 



4. My conclusions 

4.1. SPS 

Analyses of several different observables from 
independent measurements at the SPS have shown 
that a picture of production of collectively behaving 
strongly interacting matter with large volume and 
finite lifetime is now finally emerging. The 
systems produced in collisions of heaviest nuclei 
(Pb+Pb, Pb+Au) are clearly different from the 
ones produced in nuclcon-nucleon and proton- 
nucleus collisions. 

Evidence of collective behaviour of matter 
is obtained from the measured hadron yields 
and momentum spectra (NA44, NA49, WA97, 
WA98, NA50, NA52): the systematic broadening 
of transverse momentum spectra (sec. 2.1) is 
compatible with a behaviour typical for a system 
with pressure and transverse flow. The observation 
of elliptic and directed flow (sec. 2.2) also strongly 
points to this direction (NA49, WA98, NA52). The 
thermal model picture of the system reproduces 
the measured particle ratios surprisingly well 
(sec. 2.3), and, the HBT analysis of correlations in 
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SPS 


RHIC 


LHC 


y/a/A (GeV) 


17 


200 


5500 


AY 


G 


11 


17 


dN ch /dy 


400 


700-1500 


3000-8000 


r t P bpb (fm/c) 


1 


0.1 


0.005 


r QGP (We) 


1 


- 0.2 


0.1 


e(r ) (GeV/fm 3 ) 


3 


GO 


1000 


T QG p (fm/c) 


<2 


2-4 


>10 


T fo (fm/c) 


- 10 


20-30 


30-40 


Kom(T/r) 


sv Ph 


9QVp b 


400F Pb 



Table 1. Some characteristics of the heaviest systems 
produced in AA at SPS, RHIC and LHC (A ~ 200). 
From the top: cms-energy y/s, available total rapidity 
interval AY, charged particle multiplicity dN c h/dy at 
y — 0, the transit time r tr of the Lorentz-contracted 



nuclei, 



the 



formation time Tg of the QGP (lower 
limit) [112], the initial energy d ensi ty e(ro) at the 
time of formation o f the QGP |ll| , |ll2{ , lifetime of the 
QGP phase rQGP |ll2| , freeze-out time n at z = 
Jl6| , "comoving" volume V com (rf r ) along the freeze- 
out surface. Vpb denotes the rest frame size of a lead 
nucleus. The value dN c h/dy = 400 is a measured one 



Moreover, as observed by WA97, there is a clear 
specific enhancement of multistrange hadrons: the 
yield of + Q + per negative hadrons is enhanced 
by a factor ~ 10 in Pb+Pb collisions relative to 
p+p and p+Be. Microscopic hadronic rescattering 
models imply that as multistrange hadrons are 
difficult to produce due to the high mass-thresholds, 
the strangeness increase must have an origin at 
the partonic level, before hadronization. The 
agreement of the thermal models (sec. 2.3) with 
the measured hadron abundancies suggests that it 
is very likely that the system is indeed thermalized 
already before the completion of hadronization. 

Finally, the best candidate for a direct QGP 
signature so far is the anomalous suppression 
of J/^ observed in central Pb+Pb collisions by 
NA50 (sec. 2.8). It is very difficult to explain 
the additional, anomalous, suppression of J/ty 
(i.e. the suppression in addition to the expected 
suppression in nuclear matter) without invoking 
partonic degrees of freedom and deconfined matter. 
I would like to emphasize that in any possible 
alternative scenario, all the previously mentioned 
features of a collective system and the importance 
of partonic degrees of freedom will also have to be 
explained simultaneously. To my best knowledge, 
such an alternative interpretation does not exist so 
far. 

My conclusion of the results from URHIC at the 
SPS therefore is that collective behaviour of matter 
is observed but so far there is not yet a 100% certain 
proof of an observation of the QGP, although the 
strangeness enhancements and the anomalous J 
suppression strongly suggest it. 



the momentum spectra (NA44, NA49, WA98) of 
identical particles indicate that the system indeed 
is an extended one, with a large volume and hnitc 
lifetime (sec. 2.4). The pattern of the low-mass 
e + e~ enhancement, observed by CERES in S+Au 
and Pb+Au collisions (sec. 2.7), also strongly lends 
support to the picture of a collective, thermal, 
system of hadrons. The observed various pieces 
of evidence of collectivity alone do not necessarily 
imply that the system would have reached the 
QGP phase. The estimated initial energy densities 
are sufficient (table 1) but the flow studies and 
the event-by-event measurements of global variables 
(sec. 2.5) have not shown any dynamical non-trivial 
fluctuations typical for exotic events. 

However, the fact that definite additional 
anomalies are observed gives more support to the 
QGP interpretation: 

Strangeness production in A+B (sec. 2.6) is 
globally enhanced (NA35/49, WA97, NA38/50). 



4.2. SPS -> RHIC -> LHC/ ALICE 

In increasing the cms-energy of AA-collisions, the 
quantitative gains are obvious (sec. 3.2): the system 
becomes initially denser and hotter, it forms faster 
and stays collective for a longer time, i.e. its lifetime 
grows and volume increases. The longer the QGP- 
phase itself can be made, the better chances there 
are to observe direct QGP signatures. 

The higher the cms-energy is, the higher are 
also the rates for the truly hard probes of QGP 
which cannot be studied at the SPS. At RHIC 
and LHC/ ALICE, these additional QGP probes 
include T (especially its suppression) and open 
beauty production (at the LHC in particular), jets, 
high-pT direct photons, and high mass dileptons. 
Measurements of Z° may also be possible in Pb+Pb 
collisions at the CMS. At the high energies of RHIC 
and especially at the LHC/ ALICE, the hot initial 
conditions of the QGP are also expected to become 
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more reliably calculable based on perturbative 
QCD. 

According to the official plans at the time of this 
meeting, RHIC at BNL should be in operation in 
December 1999, so the first exciting data on global 
variables at high energies will be available in year 
2000. A bit further in the future, LHC/ ALICE 
at CERN will start in 2005, completing the 
extensive experimental program of ultrarelativistic 
heavy ion collisions. As the URHIC program 
at the SPS is currently teaching us, at RHIC 
and ALICE the observation of the QGP will be 
confirmed by studying carefully several independent 
observables simultaneously. After the results from 
LHC/ ALICE, a broad range of the QCD phase 
diagram has been experimentally probed, and more 
light has been shed to the theory of elementary 
particle condensed matter physics and to the early 
evolution of our Universe. 
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